arXiv:1503.00980vl [cs.AI] 3 Mar 2015 


On memetic search for the max-mean dispersion 

problem 

Xiangjing Lai and Jin-Kao Hao* 

LERIA, Universite d’Angers, 2 Bd Lavoisier, 49045 Angers, Cedex 01, France 

Submitted for publication on 19 January 2015 


Abstract 

Given a set V oi n elements and a distance matrix [dij\nxn among elements, the 
max-mean dispersion problem (MaxMeanDP) consists in selecting a subset M from 
V such that the mean dispersion (or distance) among the selected elements is maxi¬ 
mized. Being a useful model to formulate several relevant applications, MaxMeanDP 
is known to be NP-hard and thus computationally difficult. In this paper, we present 
a highly effective memetic algorithm for MaxMeanDP which relies on solution re¬ 
combination and local optimization to find high quality solutions. Computational 
experiments on the set of 160 benchmark instances with up to 1000 elements com¬ 
monly used in the literature show that the proposed algorithm improves or matches 
the published best known results for all instances in a short computing time, with 
only one exception, while achieving a high success rate of 100%. In particular, we 
improve 59 previous best results out of the 60 most challenging instances. Results 
on a set of 40 new large instances with 3000 and 5000 elements are also presented. 
The key ingredients of the proposed algorithm are investigated to shed light on how 
they affect the performance of the algorithm. 

Keywords'. Dispersion Problem; Memetic Algorithm; Tabu Search; Heuristics. 


1 Introduction 


Given a weighted complete graph G = {V,E,D), where V is the set of n 
vertices, E is the set of edges, and D represents the set of edge weights 

dij {i 7^ j), the generic equitable dispersion problem consists in selecting a 
subset M from V such that some objective function / dehned on the subgraph 
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induced by M is optimized [20]. In the related literature, a vertex n G 1^ is 
also called an element, and the edge weight dij & D is called the distance (or 
diversity) between elements i and j. 


According to the objective function to be optimized as well as the constraints 
on the cardinality of subset M, several specihc equitable dispersion problem 
can be dehned. At hrst, if the cardinality of M is hxed to a given number 
m, the related equitable dispersion problems include the following four clas¬ 
sic variants: (1) the max-sum diversity problem, also known as the maximum 
diversity problem (MDP), which is to maximize the sum of distances among 


the selected elements 


(2) the max-min diversity problem 


that aims to maximize the minimum distance among the selected elements 
[3)1 19 112 2 1123] : (3) the maximum minsum dispersion problem (MaxMinsumDP) 
that aims to maximize the minimum aggregate dispersion among the selected 
elements [3)120] : (4) the minimum differential dispersion problem (MinDiffDP) 
whose goal is to minimize the difference between the maximum and minimum 
aggregate dispersion among the selected elements to guarantee that each se¬ 
lected element has the approximately same total distance from the other se¬ 
lected elements [31I81I20] . In addition, when the cardinality of subset M is not 
hxed, i.e., the size of M is allowed to vary from 2 to n, the related equitable 
dispersion problems include the max-mean dispersion problem (MaxMeanDP) 
and the weighted MaxMeanDP 


In this study, we focus on MaxMeanDP which can be described as follows 
|20) . Given a set V oi n elements and a distance matrix {dij)nxn where dij 
represents the distance between elements i and j and can take a positive or 
negative value, the max-mean dispersion problem consists in selecting a subset 
M (|M| is not hxed) from V such that the mean dispersion among the selected 

. d - ■ 

elements, i.e., —maximized. 


MaxMeanDP can be naturally expressed as a quadratic integer program with 
binary variables x* that takes 1 if element i is selected and 0 otherwise [M], 


I.e., 


Maximize f(s) = - ^ - 

Si=l 

n 

Subject to E Xi>2 

i=l 


( 1 ) 


( 2 ) 


Xj G {0, l},i = 1,2,... ,n; (3) 

where the constraint ([2]) guarantees that at least two elements are selected. 


In addition to its theoretical signihcation as a strongly NP-hard problem |20) . 
MaxMeanDP is notable for its ability to model a variety of real-world ap¬ 
plications, such as web pages ranks na, community mining [25], and others 
mentioned in [1]. 
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Given the interest of MaxMeanDP, several useful solution approaches have 
been proposed in the literature to deal with this hard combinatorial opti¬ 
mization problem. First of all, Prokopyev et al. presented a mixed-integer 0-1 
linear programming formulation and solved small instances with up to 100 
elements with the CPLEX solver [20]. In the same work, they also introduced 
a GRASP method and made a comparison between the GRASP method and 
the MILP method to show the superiority of GRASP. Subsequently, Marti 
and Sandoya proposed a GRASP with the path relinking method (GRASP- 
PR) [TB], and the computational results show that GRASP-PR outperforms 
the previously reported methods. Very recently, Della et al. reported a hybrid 
heuristic approach based on a quadratic knapsack formulation |6], and their 
computational experiment shows that the hybrid approach is superior to the 
GRASP-PR method. In another very recent work, Garrasco et al. proposed 
a diversihed tabu search algorithm by combining a short-term tabu search 
procedure and a long-term tabu search procedure |3], and the computational 
results show that this algorithm clearly dominates the previous GRASP-PR 
method. 


In this paper, we propose the hrst population-based memetic algorithm for 
solving MaxMeanDP (called MAMMDP). The proposed algorithm combines 
a random crossover operator to generate new offspring solutions and a tabu 
search method to hnd good local optima. 


The performance of our algorithm is assessed on a set of 160 benchmark in¬ 
stances (20 < n < 1000) commonly used in the literature and a set of addi¬ 
tional 40 large-sized instances that we generate (n = 3000, 5000). For the hrst 
set of existing benchmarks, the experimental results show that the proposed 
algorithm is able to attain, in a short or very short computing time, all cur¬ 
rent best known results established by any existing algorithms, except for one 
instance. Furthermore, it can even improve the previous best known result for 
a number of these instances. The effectiveness of the proposed algorithm is 
also verihed on much larger instances of the second set with 3000 and 5000 
elements. 


The remaining part of the paper is organized as follows. In Section [H we 
describe in detail the general scheme and the ingredients of the proposed 
algorithm. In Sectional we present the computational results based on the 200 
benchmark instances and compare them with those of the existing state-of-the- 
art algorithms from the literature. In Section 0] some important ingredients of 
the proposed algorithm are analyzed and discussed. Finally, we conclude the 
paper in the last Section. 
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2 Memetic Algorithm for Max-Mean Dispersion Problem 


Memetic search is a well-known metaheuristic framework which aims to pro¬ 
vide the search with a desirable trade-off between intensification and diversi¬ 
fication through the combined use of a crossover operator (to generate new 
promising solutions) and a local optimization procedure (to locally improve 
the generated solutions) [T7IIT5] . 


2.1 General Procedure 


The proposed memetic algorithm (denoted by MAMMDP) adopts the princi¬ 
ples and guidelines of designing effective MA for discrete combinatorial prob¬ 
lems jT3]. Indeed, the overall performance of a memetic algorithm depends 
largely on the implementation of these two key components that need to be 
carefully designed according to the specific problem structure. Additionally, 
the proposed MAMMDP algorithm uses a population scheme which borrows 
from the path relinking method m to ensure a strong intensification search. 

The general procedure of our MAMMDP algorithm is shown in Algorithm [H 
where s* and respectively represent the best solution found so far and the 
worst solution in the population in terms of the objective value, and Pair Set is 
the set of solution pairs (s*, s-^), which is initially composed of all the possible 
solution pairs (s*,s-^) in the population and is dynamically updated as the 
search progresses. 

Our MAMMDP algorithm starts with an initial population P (line 4) which 
includes p different solutions, where each of them is randomly generated and 
then improved by the tabu search procedure. After the initialization of popula¬ 
tion fSection 12.3p . the algorithm enters a while loop (lines 11 to 25) to make a 
number of generations. At each generation, a solution pair (s*, s^) is randomly 
selected from PairSet and then the crossover operator (line 14) is applied to 
the selected solution pair (s*,s-^) to generate a new solution s° fSection 12.5p . 
Subsequently, s° is improved by the tabu search procedure (line 15) (Section 
12.4p . After that, a population updating rule is used to update the population 
(lines 20 to 24) fSection 12.6p . Meanwhile, the PairSet is accordingly updated 
as follows; First, the solution pair (s*,s-^) is removed from PairSet (line 13); 
Then, if an offspring solution s° replaces the worst solution s"' in the popu¬ 
lation, all the solution pairs containing are removed from PairSet and all 
the solution pairs that can be generated by combining s° with other solutions 
in the population are added into PairSet (lines 23 to 24). The while loop 
ends when PairSet becomes empty, then the population is recreated, while 
preserving the best solution (s*) found so for in the new population (lines 4 


4 




to 8), and the above while loop is repeated if the timeout limit is not reached. 

It is worth noting that compared with the traditional random selection scheme, 
the proposed MAMMDP algorithm uses the set PairSet to contain the solu¬ 
tion pairs of the population for crossover operations. This strategy ensures that 
every pair of solutions in the population is combined exactly once, favoring an 
more intensihed search. 

Algorithm 1 Memetic algorithm for Max-mean Dispersion Problem 
1: Input: The set V = {vi,V 2 , ■ ■ ■ ,Vn\ of re elements and the distance matrix D = 
[dp]nxn, the population size p, the timeout limit tout- 
2: Output: the best solution s* found 
3: repeat 

4: P = {s^,..., ^ Population_lnitialization(P) /* Section 12.31 */ 

5: if it is not in the first loop then 

6: arg mzre{/(s*) : i = 1,... ,p} 

7: P ^ P[J{s*}\{s'^} 

8: end if 

9: s* <r- arg max{/(P) : i = 1,... ,p} /* s* keeps the best solution found 

*/ 

10: PairSet ^ {(P, s-^) : 1 < f < j < p} 

11: while PairSet ^ 0 and time < tout do 

12: Randomly pick a solution pair € PairSet 

13: PairSet PairSet \ {(s*, s-^)} 

14: s° CrossoverOperator{s^, s^) /* Section 12.51 */ 

15: s° ^ TabuSearch{s°) /* Section 12.41 */ 

16: if f{s°) > f{s*) then 

17: s* ^ P 

18: end if 

19: s"' ■(— arg min{f{s^) : i = 1,... ,p} 

20: if s° dose not exist in P and f{s°) > /(s"') then 

21: P ^ PU{s°}\{s^} 

22: PairSet ^ PairSet \ {(«"', s^) : € P} 

23: PairSet PairSet U {(s°, s^) : G P} 

24: end if /* Section [2?6] */ 

25: end while 

26: until time > tout 


2.2 Search Space and Solution Representation 


Given a MaxMeanDP instance denoted by a set V of n elements as well as its 
distance matrix D = [dij]nxn, fhe search space D explored by our MAMMDP 
algorithm is composed of all possible subsets of V, i.e, D = {M : M ^ V}. 
Formally, a subset M of V can be expressed by a n-dimensional binary vector, 
(xi, X 2 ,..., Xn), where Xi takes 1 if element i belongs to M, and 0 otherwise. 
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In other words, the search space II is composed of all possible n-dimensional 
binary vectors, i.e., 

Vt = {{xi,X 2 ,... ,Xn) ■■ Xi e {0, l},l <i <n} 

Clearly, the size of the search space H is bonnded by 0(2”). 

For any candidate solntion s = (xi,X2, G II, its quality is given by the 

objective value (/(s). Formula ([T])) of the max-mean dispersion problem. 


2.3 Population Initialization 


In our memetic algorithm, the initial population of p solutions is generated 
as follows. First, we generate p random solutions, where each component xt 
(I = 1, 2,..., n) of a solution {xi,X 2 ,..., Xn) is randomly assigned a value from 
{0,1} using a uniform probability distribution. Then, the tabu search method 
(see Section [23]) is applied to each of the generated solutions to optimize them 
to a local optimum solution, and the resulting solutions are used to form the 
initial population. 


2-4 Local Optimization using Tabu Search 


Local optimization is a key component of a memetic algorithm and ensures 
generally the role of an intensihed search to locate high quality local optimum. 
In this study, we devise a tabu search (TS) method as the local optimization 
procedure which proves to be highly effective when when it is applied alone. 

Given a neighborhood structure {N{s)) and a starting solution (sq), our tabu 
search procedure iteratively replaces the incumbent solution s by a best eligible 
neighboring solution (s ) until the stopping condition is met, i.e., the best 
solution [sb) is not improved for a consecutive iterations (called the depth of 
TS). At each iteration of TS, the performed move is recorded in the tabu list to 
prevent the reverse move from being performed for the next tt iterations. Here, 
tt is called the tabu tenure and controlled by a special tabu list management 
strategy. Note that in this tabu search method a move is identihed to be 
eligible if it is not forbidden by the tabu list or it leads to a solution better 
than the best solution found so far in terms of the objective function value 
(aspiration criterion). 

The general scheme of our TS method is described in Algorithm [2l and the 
neighborhood structure employed by our TS method and the tabu list man¬ 
agement strategy are described in the following subsections. 
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Algorithm 2 TabuSearch{so, N{s),a) 

1: Input: Input solution sq. neighborhood N{s), search depth a 
2: Output: The best solution Sf, found during the tabu search process 
3: s •(— So /* s is the current solution */ 

4: Sfe ■(— s /* Sb is the best solution found so far */ 

5: d = 0 /* d counts the consecutive iterations where Sb is not updated */ 

6: repeat 

7: Choose a best eligible neighboring solution s' G A(s) 

8: s i — s' 

9: Update tabu list 

10: if /(s) > f{sb) then 

11: Sb ^— s, 

12: d = Q 

13: else 

14: d = d + l 

15: end if 

16: until d = a 
17: return Sb 


2.4-1 Move and Neighborhood 

The neighborhood Ni of our tabu search algorithm is dehned by the one-flip 
move operator which consists of changing the value of a single variable x* 
to its complementary vale 1 — x*. As such, given a solution s, the one-flip 
neighborhood A^i(s) of s is composed of all possible solutions that can be 
obtained by applying the one-flip move to s. The size of the neighborhood 
Ni{s) is thus bounded by 0{n), where n is the number of components in s. 


2 . 4.2 Fast Neighborhood Evaluation Teehnique 

In our TS method, we employ a fast neighborhood evaluation technique to 
examine the neighborhood Ni. For this purpose, we maintain a n-dimensional 
vector W = (pi,P2, • • • ,Pn) to effectively calculate the move value (i.e., the 
change of objective value) of each possible move applicable to the current so¬ 
lution, where the entry pi represents the sum of distances between the element 
i and the selected elements for the current solution, i.e., Pi = 
where M is the set of selected elements. 

Given the current solution s, if an one-flip move is performed, i.e., a variable 
Xj is flipped as Xj (1 — Xj), then the move value A* can be rapidly computed 
as follows: 
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Ai = < 


Pi 


for Xi = 0 ; 


( 4 ) 


-f{s) 


|M| + 1 

fis) 


+ 


\M\ 


|M| + 1’ 

Pi 

\M\ - 1’ 


for Xi = 1 ; 


( 5 ) 


where /(s) is the objective value of the current solution s and \M\ is the 
number of selected elements in s. Subsequently, the vector W is accordingly 
updated as: 

r pj + dij, for Xi = 0J 7^ i-, ( 6 ) 

Pj = \ Pj- dij, for Xi = l,j ^ i] (7) 

[pj, ioTj = i; (8) 

Note that the vector W is initialized at the beginning of each call of TS with 

the complexity of O(n^), and is updated in 0(n) after each move. 


2.4-3 Tabu List Management Strategy 

In our TS procedure, we use a tabu list management strategy to dynamically 
tune the tabu tenure tt, which is adapted according to a technique proposed 
in where the tabu tenure is given by a periodic step function. If the current 
iteration is y, then the tabu tenure of a move is denoted by tt{y). 

Precisely, our tabu tenure function is defined, for each period, by a sequence of 
values (ui, 02, • • • , Oq+i) and a sequence of interval margins (?/i, ?/2, •'' ) Vq+i) 
such that for each y in — 1 ], tt{y) = Oj + rand{2), where rand{2) de¬ 

notes a random integer between 0 to 2. Here, q is hxed to 15, (a)i=i,... ,15 = 
2 , 1 , 4, 1 , 2 , 1 , 8 , 1 , 2 , 1 , 4, 1 , 2 , 1 ), where T^ax is a parameter which rep¬ 
resents the maximum tabu tenure. Finally, the interval margins are defined by 
yi = 1, yi+i = yi + 5ai {i < 15). 


Thus, this function varies periodically and for each period, 15 tabu tenures 
are used dynamically, each being kept for a number of consecutive iterations. 
In principle, this function helps the tabu search procedure reach a desirable 
tradeoff between intensification and diversification during its search. 


2.5 Crossover Operator 


In a memetic algorithm, the crossover operator is another essential ingredient 
whose main goal is to bring the search process to new promising search regions 
to diversify the search. In this work, we investigate two crossover operators 
for MaxMeanDP; the first one is the standard uniform crossover (denoted by 
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Algorithm 3 The Uniform Crossover Operator for MaxMeanDP 
1: Input: Two parent solutions = {x\,x\, ..., x\) and = (xf, , x^). 

2: Output: Offspring solution s° = (x°,X 2 ,.. • ,x") 

3: for i = 1 to n do 

4: r ^ rand[Q, 1) /* rand[0,l) denotes a random number between 0 and 1 */ 

5: if r < 0.5 then 

6 : x° ^ x] 

7: else 

8: x° ^ X? 

9: end if 

10 : end for 

11 : return = (x°, X 2 ,..., x°) 


Algorithm 4 The Greedy Crossover Operator for MaxMeanDP 
1: Input: Two parent solutions and and their subsets of selected elements are 
respectively denoted by Mi and M 2 

2: Output: Offspring solution s° whose subset of selected elements is denoted by Mq 
3: m ■<— (|Mi| + |M2|)/2 /* Determine approximately the size of the set Mq */ 

/*Generate a partial solution by preserving the common elements of Mi and M 2 */ 

4: Mo ^ Ml n M 2 

5: M[ ^ Ml \ Mo, M 2 ^ M 2 \ Mo 

6 : while |Mo| < m do 

7: if M{ 7 ^ 0 then 

8: Xo argmax{IS.f{v) : v € M;^} 

/* A.f{v) represents the move value of adding element v to Mo */ 

9: Mo ^ Mo U {ro}, M[ ^ M[\ {vo} 

10 : end if 

11 : if Mg / 0 then 

12: Uo ^ argmax{Af{u) : u € Mg} 

13: Mo ^ Mo U {no}, M'^ ^ M'^\ {uo} 

14: end if 

15: end while 

16: return s°(i.e.,Mo) 

UC) operator, and the other is a specific greedy crossover (denoted by GC) 
operator. 

UC is very simple and is described in Algorithm [3l Given two parent solntions 
= (xj, Xg,..., and = (xf, Xg,..., x^), the valne of each component 
x° {i = 1,2,... ,77.) of the offspring solntion s° is randomly chosen from the 
set {x}, x^} with the same probability of 0.5. In spite of its simplicity, UC has 
shown to be qnite robnst and effective in many settings. 

The dedicated GC operator is shown in Algorithm HI Given two parent soln¬ 
tions and whose snbsets of selected elements are respectively represented 
by Ml and M 2 , we first estimate the size of the set Mo of selected elements 
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for the offspring solution s° such that the Hamming distances of s° to as 
well as are approximately the same, i.e., m ■(— Then a partial 

solution is generated by preserving the common elements of Mi and M 2 , i.e., 
Mo •(— Ml n M 2 . After that, we complete the offspring solution s° in a step by 
step way by choosing in turn a best element from Mi \ Mo or M 2 \ Mo and 
adding it to Mo- More specifically, we first choose an element {vo) that yields 
the largest move value from Mi \ Mo and add it to Mo, then we choose another 
element {uo) that yields the largest move value from M 2 \ Mo and add it to 
Mo- We repeat these two steps until the size of Mo reaches m. 

Intuitively, UC is more disruptive than GC and thus is suitable for the pur¬ 
pose of diversifying the search. On the other hand, GC is more conservative 
and computationally more expensive. Based on the computational outcomes 
presented in Section we adopt in this study the UC operator as the main 
crossover operator of our memetic algorithm, while using GC as a reference 
operator for our analysis on crossovers. 


2.6 Population Updating Rule 


When a new offspring solution is generated by the crossover operator, it is first 
improved the tabu search procedure and then used to update the population 
according to the following rule. If the offspring solution is distinct from any 
existing solution in the population and is better than the worst solution in the 
population in terms of objective value, then the offspring solution replaces the 
worst solution of the population. Otherwise, the population is kept unchanged. 


3 Experimental Results and Comparisons 


In this section, we run extensive computational experiments to assess the per¬ 
formance of our memetic algorithm based on a large number of MaxMeanDP 
benchmark instances. 


3.1 Benchmark Instances 


Our computational experiments are carried out on two types of instances, 
namely Type I and Type II. The distances of Type I instances are randomly 
generated in the interval [—10,10] with a uniform probability distribution, 
while the distances of Type II instances are generated from [—10, —5] U [5,10] 
with the same probability distribution. 
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Additionally, the set of benchmark instances used in our experiments is com¬ 
posed of two subsets. The hrst subset consists of 80 Type I instances and 80 
Type II instances with the number of elements n ranging from 20 to 1000. 
These 160 instances were extensively adopted by the previous studies 
and are available online at http: //www. optsicom. es. The second subset con¬ 
sists of 20 Type I and 20 Type II large instances with n = 3000 or 5000. The 
source code of the generator used to obtain these 40 large instances will be 
availabl^. 


3.2 Parameter Settings and Experimental Protocol 


Our memetic algorithm relies on only three parameters: the population size 
p, the depth of tabu search a and the maximum tabu tenure Tmax- For p 
and a, we follow j21] and set p = 10, a = 50000 while setting T^ax = 120 
empirically. This parameter setting is used for all the experiments reported in 
the paper. Even if hne-tuning these parameters would lead to better results, 
as we show below, our algorithm with this hxed setting is able to attain a high 
performance with respect to the state of the art results. 


Our memetic algorithm is programmed in C++ and compiled using g++ com¬ 
piler with the ’-02’ fla^ZI- All experiments are carried out on a computer with 
an Intel Xeon E5440 processor (2.83 GHz CPU and 2Gb RAM), running the 
Linux operating system. Following the DIMACS machine benchmark proce¬ 
dure, our machine requires respectively 0.23, 1.42, and 5.42 seconds for the 
graphs r300.5, r400.5, r500.5. 


Given the stochastic nature of our algorithm, we solve each tested problem 
instance 20 times, where the stopping condition is given by a cutoff time limit 
which depends on the size of the instances. Specihcally, the cutoff limit tout is 
set to be 10 seconds for n < 150, 100 seconds for n G [500,1000], 1000 seconds 
for n = 3000, and 2000 seconds for n = 5000. As we discuss in Section 13.31 
these time limits are signihcantly shorter than those used by the reference 
algorithms of the literature. 
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Computational results of the proposed MAMMDP algorithm on the set of 60 rep¬ 
resentative instances with 500 < n < 1000. Each instance is independently solved 
20 times, and improved results are indicated in bold compared to the previous best 
known results /p^e of the literature reported in mm- 


Instance 

n 

fpre 

lUFTTSl 


Memetic Algori 

thm 


Sbest 

lavg 

yH 


MUPll 500 

SCD 

81.28 

81.277044 

81.277044 

20/20 

CTBS 

MDPI2 500 

500 

77.60 

78.610216 

78.610216 

20/20 

1.43 

MDPI3 500 

500 

75.65 

76.300787 

76.300787 

20/20 

2.71 

MDPI4 500 

500 

82.28 

82.332081 

82.332081 

20/20 

0.95 

MDPI5 500 

500 

80.01 

80.354029 

80.354029 

20/20 

2.80 

MDPI6 500 

500 

81.12 

81.248553 

81.248553 

20/20 

0.78 

MDPI7 500 

500 

78.09 

78.164511 

78.164511 

20/20 

0.92 

MDPI8 500 

500 

79.01 

79.139881 

79.139881 

20/20 

1.27 

MDPI9 500 

500 

77.15 

77.421000 

77.421000 

20/20 

2.37 

MDPIIO 500 

500 

81.24 

81.309871 

81.309871 

20/20 

0.91 

MDPIIl 500 

500 

109.33 

109.610136 

109.610136 

20/20 

0.75 

MDPII2 500 

500 

105.06 

105.717536 

105.717536 

20/20 

0.88 

MDPII3 500 

500 

107.64 

107.821739 

107.821739 

20/20 

0.89 

MDPII4 500 

500 

105.69 

106.100071 

106.100071 

20/20 

0.56 

MDPII5 500 

500 

106.59 

106.857162 

106.857162 

20/20 

0.99 

MDPII6 500 

500 

106.17 

106.297958 

106.297958 

20/20 

0.98 

MDPII7 500 

500 

106.92 

107.149379 

107.149379 

20/20 

0.88 

MDPII8 500 

500 

103.49 

103.779195 

103.779195 

20/20 

0.59 

MDPII9 500 

500 

106.20 

106.619793 

106.619793 

20/20 

1.11 

MDPIIIO 500 

500 

103.79 

104.651507 

104.651507 

20/20 

1.01 

MDPIl 750 

750 

95.86 

96.644236 

96.644236 

20/20 

7.98 

MDPI2 750 

750 

97.42 

97.564880 

97.564880 

20/20 

3.82 

MDPI3 750 

750 

96.97 

97.798864 

97.798864 

20/20 

1.81 

MDPI4 750 

750 

95.21 

96.041364 

96.041364 

20/20 

4.38 

MDPI5 750 

750 

96.65 

96.740448 

96.740448 

20/20 

1.21 

MDPI6 750 

750 

99.25 

99.861250 

99.861250 

20/20 

5.55 

MDPI7 750 

750 

96.26 

96.545413 

96.545413 

20/20 

1.01 

MDPI8 750 

750 

96.46 

96.726976 

96.726976 

20/20 

1.73 

MDPI9 750 

750 

96.78 

98.058377 

98.058377 

20/20 

2.18 

MDPIIO 750 

750 

99.85 

100.064185 

100.064185 

20/20 

3.42 

MDPIIl 750 

750 

127.69 

128.863707 

128.863707 

20/20 

5.66 

MDPII2 750 

750 

130.79 

130.954426 

130.954426 

20/20 

2.31 

MDPII3 750 

750 

129.40 

129.782453 

129.782453 

20/20 

11.64 

MDPII4 750 

750 

125.68 

126.582271 

126.582271 

20/20 

1.48 

MDPII5 750 

750 

128.13 

129.122878 

129.122878 

20/20 

1.32 

MDPII6 750 

750 

128.55 

129.025215 

129.025215 

20/20 

7.98 

MDPII7 750 

750 

124.91 

125.646682 

125.646682 

20/20 

3.38 

MDPII8 750 

750 

130.66 

130.940548 

130.940548 

20/20 

1.91 

MDPII9 750 

750 

128.89 

128.889908 

128.889908 

20/20 

1.30 

MDPIIIO 750 

750 

132.99 

133.265300 

133.265300 

20/20 

1.81 

MDPIl 1000 

1000 

118.76 

119.174112 

119.174112 

20/20 

8.25 

MDPI2 1000 

1000 

113.22 

113.524795 

113.524795 

20/20 

3.52 

MDPI3 1000 

1000 

114.51 

115.138638 

115.138638 

20/20 

2.32 

MDPI4 1000 

1000 

110.53 

111.127437 

111.127437 

20/20 

5.72 

MDPI5 1000 

1000 

111.24 

112.723188 

112.723188 

20/20 

1.61 

MDPI6 1000 

1000 

112.08 

113.198718 

113.198718 

20/20 

7.72 

MDPI7 1000 

1000 

110.94 

111.555536 

111.555536 

20/20 

1.88 

MDPI8 1000 

1000 

110.29 

111.263194 

111.263194 

20/20 

3.55 

MDPI9 1000 

1000 

115.78 

115.958833 

115.958833 

20/20 

2.38 

MDPIIO 1000 

1000 

114.29 

114.731644 

114.731644 

20/20 

2.16 

MDPIIl 1000 

1000 

145.46 

147.936175 

147.936175 

20/20 

1.60 

MDPII2 1000 

1000 

150.49 

151.380035 

151.380035 

20/20 

1.78 

MDPII3 1000 

1000 

149.36 

150.788178 

150.788178 

20/20 

4.92 

MDPII4 1000 

1000 

147.91 

149.178006 

149.178006 

20/20 

3.80 

MDPII5 1000 

1000 

150.23 

151.520847 

151.520847 

20/20 

3.28 

MDPII6 1000 

1000 

147.29 

148.343378 

148.343378 

20/20 

3.22 

MDPII7 1000 

1000 

148.41 

148.742375 

148.742375 

20/20 

6.30 

MDPII8 1000 

1000 

145.87 

147.826804 

147.826804 

20/20 

13.52 

MDPII9 1000 

1000 

145.67 

147.078145 

147.078145 

20/20 

1.63 

MDPIIIO 1000 

1000 

148.40 

150.046137 

150.046137 

20/20 

2.13 

^Better 



50 




T^Equal 



1 

1 



7 ^ Worse 



0 

0 




3.3 Computational Results and Comparisons on Small and Medium Sized 
Instances 

Our first experiment aims to evaluate the performance of our MAMMDP 
algorithm on the set of 160 popular instances with up to 1000 elements. The 


^ The source code of generating these instances will be available from our website. 
^ Our best results and the source code of our algorithm will be made available 
online. 

^ dmclique, ftp://dimacs.rutgers.edu/pub/dsj/clique, the benchmark procedure is 
complied by gcc compiler with the ’-02’ flag 
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computational results of MAMMDP on the 60 medium sized instances are 
summarized in Table [H whereas the results of the 100 small instances with 
n < 150 are given in the Appendix [6] with the same computational statistics. 

The hrst two columns of the table give respectively the name and size of 
instances. Column 3 indicates the best objective values {fpre) of the literature 
which are compiled from the best results yielded by three recent and best 
performing algorithms, namely GRASP-PR jTB], a hybrid heuristic approach 
[B], and a diversihed tabu search method |3] (from http: //www. optsicom. es). 
Note that the previous best known results {fpre) are given with two decimal 
in the literature. In [TlllG) . the cutoff time limits are set to 90, 600, and 1800 
seconds for instances with the size 500, 750, and 1000, respectively, and in 
[B] the cutoff time limits are respectively set to 60 and 600 seconds for the 
instances with the size 150 and 500. The GRASP-PR method was performed 
on a computer with an Intel Gore Solo 1.4 GHz GPU with 3 GB RAM |16| : the 
hybrid heuristic approach was run on a computer with an Intel Gore 15-3550 
3.30GHz GPU with 4GB RAM jB] and the diversihed tabu search method was 
run on a computer with an Intel Gore 2 Quad GPU and 6 GB RAM |3]. 

Our results are reported in columns 4 to 7, including the best objective value 
{fbest) yielded over 20 independent runs, the average objective value {favg), the 
success rate {SR) to achieve fbest, and the average computing time in seconds 
{t{s)) to achieve fbest- The last three rows Better, Equal, Worse of the table 
respectively show the number of instances for which our result is better, equal 
to and worse than fpre- The improved results are indicated in bold compared 

t*^ fpre- 

First, one observes from Table [T] that our MAMMDP algorithm improves 
the previous best known result for all instances except for one instance for 
which our result matches the previous best known result. Therefore, these 
results clearly indicate the superiority of the proposed MAMMDP algorithm 
compared to the previous MaxMeanDP algorithms. Second, when examining 
the success rate of the algorithm, one can hnd that the MAMMDP algorithm 
achieves a success rate of 100% for all tested instances, which means a good 
robustness of our MAMMDP algorithm. Third, in terms of average computing 
time, one observes that for all instances our MAMMDP algorithm obtains its 
best result with an average time of less than 14 seconds, which are much 
shorter than those of the previous algorithms in the literature. 


3-4 Computational Results on Large-Scale Instances 


In order to further assess the performance of the proposed MAMMDP algo¬ 
rithm on the large-scale instances, we run independently MAMMDP 20 times 
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Table 2 

Computational results of the proposed memetic algorithm on the set of 40 large 
instances with n = 3000, 5000. Each instance is independently solved 20 times. 


Memetic Algorithm 


Instance 

n 

fbest 

iavq 

Sr 

t{s) 

MDPIl 3000 

3000 

189.048965 

189.048965 

20/20 

88.36 

MDPI2 3000 

3000 

187.387292 

187.387292 

20/20 

60.71 

MDPI3 3000 

3000 

185.666806 

185.655084 

13/20 

352.85 

MDPI4 3000 

3000 

186.163727 

186.153631 

16/20 

300.37 

MDPI5 3000 

3000 

187.545515 

187.545515 

20/20 

61.29 

MDPI6 3000 

3000 

189.431257 

189.431257 

20/20 

51.99 

MDPI7 3000 

3000 

188.242583 

188.242583 

20/20 

86.57 

MDPI8 3000 

3000 

186.796814 

186.796814 

20/20 

48.04 

MDPI9 3000 

3000 

188.231264 

188.231264 

20/20 

151.78 

MDPIIO 3000 

3000 

185.682511 

185.623778 

10/20 

228.72 

MDPIIl 3000 

3000 

252.320433 

252.320433 

20/20 

59.70 

MDPII2 3000 

3000 

250.062137 

250.062137 

20/20 

220.10 

MDPII3 3000 

3000 

251.906270 

251.906270 

20/20 

146.32 

MDPII4 3000 

3000 

253.941007 

253.940596 

19/20 

370.76 

MDPII5 3000 

3000 

253.260423 

253.260350 

17/20 

374.00 

MDPII6 3000 

3000 

250.677750 

250.677750 

20/20 

55.35 

MDPII7 3000 

3000 

251.134413 

251.134413 

20/20 

74.72 

MDPII8 3000 

3000 

252.999648 

252.999648 

20/20 

79.82 

MDPII9 3000 

3000 

252.425770 

252.425770 

20/20 

90.27 

MDPIIIO 3000 

3000 

252.396590 

252.396590 

20/20 

13.18 

MDPIl 5000 

5000 

240.162535 

240.102875 

7/20 

312.13 

MDPI2 5000 

5000 

241.827401 

241.792978 

6/20 

1244.36 

MDPI3 5000 

5000 

240.890819 

240.888162 

19/20 

810.48 

MDPI4 5000 

5000 

240.997186 

240.976789 

6/20 

653.64 

MDPI5 5000 

5000 

242.480129 

242.475885 

19/20 

735.16 

MDPI6 5000 

5000 

240.322850 

240.306326 

8/20 

976.02 

MDPI7 5000 

5000 

242.814943 

242.774982 

5/20 

259.50 

MDPI8 5000 

5000 

241.194990 

241.161763 

8/20 

1148.60 

MDPI9 5000 

5000 

239.760560 

239.667613 

4/20 

1219.71 

MDPIIO 5000 

5000 

243.473734 

243.373015 

4/20 

457.28 

MDPIIl 5000 

5000 

322.235897 

322.181291 

5/20 

1519.05 

MDPII2 5000 

5000 

327.301910 

327.006342 

5/20 

1103.13 

MDPII3 5000 

5000 

324.813456 

324.801590 

10/20 

955.81 

MDPII4 5000 

5000 

322.237586 

322.197276 

5/20 

664.10 

MDPII5 5000 

5000 

322.491211 

322.380726 

7/20 

1014.90 

MDPII6 5000 

5000 

322.950488 

322.703887 

4/20 

352.88 

MDPII7 5000 

5000 

322.850438 

322.793125 

10/20 

714.31 

MDPII8 5000 

5000 

323.112120 

323.053268 

11/20 

879.48 

MDPII9 5000 

5000 

323.543775 

323.339842 

7/20 

569.73 

MDPIIIO 5000 

5000 

324.519908 

324.414458 

15/20 

752.95 


to solve each instance of the second set of benchmarks with n > 3000, and 
report the computational statistics in Table [2] with the same information as 
in Table [1] 


Table [2] discloses that for the instances with 3000 elements, our MAMMDP 
algorithm reaches a success rate of at least 10/20, which is an interesting 
indicator as to its good performance for these instances. However, for the still 
larger instances with n = 5000, the success rate of the algorithm significantly 
varies between 4/20 and 19/20, which means that these large instances are 
clearly more difficult. Nevertheless, one observes that the difference between 
the best and average objective values is very small for all instances. These 
results can be served as reference lower bounds for future comparisons of new 
MaxMeanDP algorithms. 


4 Analysis and Discussions 


In this section, we study some essential ingredients of the proposed algorithm 
to understand their impacts on the performance of the proposed algorithm. 
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4-1 Effectiveness of the Tabu Search Procedure 


Table 3 

Computational results of the underlying tabu search procedure of the proposed algo¬ 
rithm on the set of 30 representative instances with 500 < n < 1000. Each instance is 
independently solved 20 times, and improved results are indicated in bold compared 
to the previous best known objective values. 


Instance 

n 

/p.e 

/* 

fbest 

Tabu Search 
favg 

SR 

t(s) 

MDPIl 500 

500 

81.28 

81.277044 

81.277044 

81.277044 

20/20 

0.54 

MDPI2 500 

500 

77.60 

78.610216 

78.610216 

78.586017 

8/20 

0.68 

MDPI3 500 

500 

75.65 

76.300787 

76.300787 

76.252735 

7/20 

0.80 

MDPI4 500 

500 

82.28 

82.332081 

82.332081 

82.325976 

15/20 

0.65 

MDPI5 500 

500 

80.01 

80.354029 

80.354029 

80.349285 

5/20 

0.78 

MDPIIl 500 

500 

109.33 

109.610136 

109.610136 

109.473283 

15/20 

0.56 

MDPII2 500 

500 

105.06 

105.717536 

105.717536 

105.678373 

14/20 

0.70 

MDPII3 500 

500 

107.64 

107.821739 

107.821739 

107.808748 

18/20 

0.77 

MDPII4 500 

500 

105.69 

106.100071 

106.100071 

106.08738 

18/20 

0.54 

MDPII5 500 

500 

106.59 

106.857162 

106.857162 

106.834002 

14/20 

0.64 

MDPIl 750 

750 

95.86 

96.644236 

96.644236 

96.472938 

5/20 

1.24 

MDPI2 750 

750 

97.42 

97.564880 

97.564880 

97.557159 

9/20 

1.10 

MDPI3 750 

750 

96.97 

97.798864 

97.798864 

97.780265 

15/20 

1.05 

MDPI4 750 

750 

95.21 

96.041364 

96.041364 

95.880611 

2/20 

1.17 

MDPI5 750 

750 

96.65 

96.740448 

96.740448 

96.731085 

18/20 

0.99 

MDPIIl 750 

750 

127.69 

128.863707 

128.863707 

128.401188 

4/20 

1.12 

MDPII2 750 

750 

130.79 

130.954426 

130.954426 

130.920628 

7/20 

1.01 

MDPII3 750 

750 

129.40 

129.782453 

129.782453 

129.557185 

1/20 

1.03 

MDPII4 750 

750 

125.68 

126.582271 

126.582271 

126.441529 

11/20 

0.96 

MDPII5 750 

750 

128.13 

129.122878 

129.122878 

129.032812 

15/20 

1.11 

MDPIl 1000 

1000 

118.76 

119.174112 

119.174112 

118.836369 

5/20 

1.63 

MDPI2 1000 

1000 

113.22 

113.524795 

113.524795 

113.363683 

7/20 

1.74 

MDPI3 1000 

1000 

114.51 

115.138638 

115.138638 

115.106057 

16/20 

1.33 

MDPI4 1000 

1000 

110.53 

111.127437 

111.127437 

110.924380 

4/20 

1.42 

MDPI5 1000 

1000 

111.24 

112.723188 

112.723188 

112.723188 

20/20 

1.18 

MDPIIl 1000 

1000 

145.46 

147.936175 

147.936175 

147.919913 

17/20 

1.44 

MDPII2 1000 

1000 

150.49 

151.380035 

151.380035 

151.344939 

18/20 

1.54 

MDPII3 1000 

1000 

149.36 

150.788178 

150.788178 

150.475002 

7/20 

1.82 

MDPII4 1000 

1000 

147.91 

149.178006 

149.178006 

149.093037 

11/20 

1.59 

MDPII5 1000 

1000 

150.23 

151.520847 

151.520847 

151.495094 

10/20 

1.57 

T^Better 




29 

29 



:^Equal 




1 

1 



Worse 




0 

0 




First, to show the effectiveness of the underlying tabu search procedure of 
the proposed algorithm, we carry out an additional experiment on 30 rep¬ 
resentative instances, where each instance is independently solved 20 times 
by the tabu search procedure with randomly generated initial solutions. The 
computational results are summarized in Table [3l where f* and t[s) represent 
respectively the previous best known result and the average computing time 
over 20 runs, other symbols are the same as those in Tabled] Note that the im¬ 
proved results compared to the previous best known values /p^e are indicated 
in bold. 


Table [3] discloses that our tabu search procedure alone is able to attain the 
previous best known result for the tested instances. Moreover, for all tested 
instances, the obtained average objective value {favg) is still better than the 
previous best known result {fpre) on 29 out of 30 instances. These results 
indicate that our tabu search method is quite effective with respect to the 
existing algorithms designed for MaxMeanDP. In terms of computational time, 
one observes that one run of the tabu search procedure takes on average less 
than 2 seconds for instances with n < 1000, much shorter than those required 
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by the state of the art algorithms mEG\- To sum, this experiment shows 
that our tabu search procedure is highly effective compared with the existing 
state-of-the-art algorithms in the literature. 

It should be mentioned that compared to the previous methods for MaxMe- 
auDP, such as those in |1], the success of our tabu search method may be 
attributed to the combined use of the neighborhood A^i, the fast neighbor¬ 
hood evaluation technique and its dynamic tabu list management strategy. 


4-2 Influence of the Crossover Operator 


Table 4 

Comparison between the UC and GC crossover operators on the set of 40 large 
instances with n = 3000 or 5000. Each instance is independently solved 20 times 
by the memetic algorithms with the UC and GC operators respectively, and better 
results between these two memetic algorithms are indicated in bold . 


Uniform Crossover Greedy Crossover 


Instance 


fbest 

favq 

SR 

t{s) 

fbest 

favg 

SR 

t(s) 

MDPIl 

3000 

189.048965 

189.048965 

20/20 

88.36 

189.048965 

189.048965 

20/20 

82.49 

MDPI2 

3000 

187.387292 

187.387292 

20/20 

60.71 

187.387292 

187.387292 

20/20 

71.95 

MDPI3 

3000 

185.666806 

185.655084 

13/20 

352.85 

185.666806 

185.656214 

14/20 

430.55 

MDPI4 

3000 

186.163727 

186.153631 

16/20 

300.37 

186.163727 

186.163727 

20/20 

248.47 

MDPI5 

3000 

187.545515 

187.545515 

20/20 

61.29 

187.545515 

187.545515 

20/20 

87.26 

MDPI6 

3000 

189.431257 

189.431257 

20/20 

51.99 

189.431257 

189.431257 

20/20 

51.15 

MDPI7 

3000 

188.242583 

188.242583 

20/20 

86.57 

188.242583 

188.242583 

20/20 

103.65 

MDPI8 

3000 

186.796814 

186.796814 

20/20 

48.04 

186.796814 

186.796814 

20/20 

68.57 

MDPI9 

3000 

188.231264 

188.231264 

20/20 

151.78 

188.231264 

188.231264 

20/20 

69.86 

MDPIIO 

3000 

185.682511 

185.623778 

10/20 

228.72 

185.682511 

185.623041 

14/20 

386.07 

MDPIIl 

3000 

252.320433 

252.320433 

20/20 

59.70 

252.320433 

252.320433 

20/20 

78.92 

MDPII2 

3000 

250.062137 

250.062137 

20/20 

220.10 

250.062137 

250.062137 

20/20 

219.11 

MDPII3 

3000 

251.906270 

251.906270 

20/20 

146.32 

251.906270 

251.906270 

20/20 

139.13 

MDPII4 

3000 

253.941007 

253.940596 

19/20 

370.76 

253.941007 

253.938635 

19/20 

243.89 

MDPII5 

3000 

253.260423 

253.260350 

17/20 

374.00 

253.260423 

253.260423 

20/20 

301.14 

MDPII6 

3000 

250.677750 

250.677750 

20/20 

55.35 

250.677750 

250.677750 

20/20 

63.01 

MDPII7 

3000 

251.134413 

251.134413 

20/20 

74.72 

251.134413 

251.134413 

20/20 

75.89 

MDPII8 

3000 

252.999648 

252.999648 

20/20 

79.82 

252.999648 

252.999648 

20/20 

62.52 

MDPII9 

3000 

252.425770 

252.425770 

20/20 

90.27 

252.425770 

252.425770 

20/20 

96.89 

MDPIIIO 3000 

252.396590 

252.396590 

20/20 

13.18 

252.396590 

252.396590 

20/20 

18.34 

MDPIl 

5000 

240.162535 

240.102875 

7/20 

312.13 

240.162535 

240.110086 

9/20 

104.32 

MDPI2 

5000 

241.827401 

241.792978 

6/20 

1244.36 

241.827401 

241.750738 

3/20 

873.37 

MDPI3 

5000 

240.890819 

240.888162 

19/20 

810.48 

240.890819 

240.875111 

14/20 

681.51 

MDPI4 

5000 

240.997186 

240.976789 

6/20 

653.64 

240.997186 

240.987885 

10/20 

1062.58 

MDPI5 

5000 

242.480129 

242.475885 

19/20 

735.16 

242.480129 

242.475876 

18/20 

790.37 

MDPI6 

5000 

240.322850 

240.306326 

8/20 

976.02 

240.376038 

240.296124 

4/20 

564.79 

MDPI7 

5000 

242.814943 

242.774982 

5/20 

259.50 

242.820139 

242.782093 

5/20 

261.98 

MDPI8 

5000 

241.194990 

241.161763 

8/20 

1148.60 

241.194990 

241.159647 

10/20 

1219.58 

MDPI9 

5000 

239.760560 

239.667613 

4/20 

1219.71 

239.760560 

239.569628 

4/20 

756.59 

MDPIIO 

5000 

243.473734 

243.373015 

4/20 

457.28 

243.385487 

243.356904 

9/20 

1057.73 

MDPIIl 

5000 

322.235897 

322.181291 

5/20 

1519.05 

322.235897 

322.172704 

5/20 

541.47 

MDPII2 

5000 

327.301910 

327.006342 

5/20 

1103.13 

327.301910 

326.953531 

7/20 

662.58 

MDPII3 

5000 

324.813456 

324.801590 

10/20 

955.81 

324.813456 

324.790747 

4/20 

629.53 

MDPII4 

5000 

322.237586 

322.197276 

5/20 

664.10 

322.237586 

322.205163 

7/20 

535.54 

MDPII5 

5000 

322.491211 

322.380726 

7/20 

1014.90 

322.491211 

322.403539 

8/20 

1095.86 

MDPII6 

5000 

322.950488 

322.703887 

4/20 

352.88 

322.950488 

322.692191 

6/20 

457.37 

MDPII7 

5000 

322.850438 

322.793125 

10/20 

714.31 

322.850438 

322.774304 

12/20 

631.95 

MDPII8 

5000 

323.112120 

323.053268 

11/20 

879.48 

323.112120 

323.007841 

13/20 

764.65 

MDPIIO 

5000 

323.543775 

323.339842 

7/20 

569.73 

323.543775 

323.273630 

5/20 

305.89 

MDPIIIO 5000 

324.519908 

324.414458 

15/20 

752.95 

324.519908 

324.500667 

19/20 

659.21 

T^Better 


1 

16 



2 

9 



T^Equal 


37 

15 



37 

15 



9 ^ Worse 


2 

9 



1 

16 



p-value 


5.637e-l 

1.615e-l 








In this section, we show a study about the influence of the crossover opera¬ 
tor on the performance of the proposed algorithm by comparing the uniform 
crossover (UC) operator and the greedy crossover (GC) operator. The exper¬ 
iment is carried out on a set of 40 large instance with n > 3000, where the 
memetic algorithms with UC and GC are respectively performed 20 times 
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for each of the tested instances. The computational results of the both al¬ 
gorithms are summarized in Table IU including the best {fbest) and average 
ifavg) objective values obtained over 20 runs, the success rate {SR) and av¬ 
erage computing time in seconds (t(s)) to achieve the associated fbest, where 
better results between these two algorithms are indicated in bold. The rows, 
Better, Equal, Worse denote respectively the number of instances for which 
an algorithm yields better, equal, worse results compared to the other algo¬ 
rithm. Finally, to verify whether there exists a significant difference between 
the two crossover operators in terms of the best and average objective values, 
the p-values from the non-parametric Friedman test are reported in the last 
row. 

One observes from Table 0] that the UC and GC operators are comparable in 
the overall performance of the algorithm in terms of best and average results, 
which is confirmed by the Friedman test (p-values > 0.05). First, in terms of 
the best objective value, UC and GC yield respectively better result on 1 and 
2 instances compared to another operator. Second, in terms of the average 
objective value, it can be found the UC operator produces a better and worse 
result respectively for 16 and 9 instances. As to the success rate and computing 
time, these two crossover operators also achieve similar performances. This 
experiment demonstrates that there is no dominance of one crossover over 
the other. Instead, they are complementary to solve different instances. One 
interesting future study would be to investigate the ways of using jointly these 
two operators with the search algorithm. 


4-3 Improvement of Memetic Algorithm Over the Tabu Search Procedure 


As shown in Section 14.11 our tabu search procedure is very competitive com¬ 
pared to the existing algorithms in the literature. So it is interesting to know 
whether our MAMMDP algorithm has a significant improvement over its un¬ 
derlying local optimization (the tabu search) procedure. To compare the per¬ 
formances of the MAMMDP algorithm and its underlying tabu search pro¬ 
cedure, the multi-start tabu search (MTS) and MAMMDP algorithms are 
respectively performed 20 times for each of the 40 representative instances 
with n = 3000 or 5000 under the same cutoff time limits given in Section 13.21 
Notice that for MTS, the tabu search procedure is run in a multi-start way 
with a randomly generated initial solution for each re-start until the timeout 
limit is reached, the tabu search procedure being re-started once the depth of 
tabu search a (which is set to 5 x 10^) is reached. The computational results 
of both algorithms are respectively summarized in Table [5] which is composed 
of two parts, where the symbols are the same with those in Table 01 

Table 0] discloses that for the 20 instances with n = 3000 the MAMMDP 
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Table 5 

Comparison between the multi-start tabu search method (MTS) and the proposed 
memetic algorithm on the set of 40 large instances with n > 3000. Each instance 
is independently solved 20 times by both algorithms respectively, and better results 
between two algorithms are indicated in bold. 


Memetic Algorithm MTS 


Instance 


fbest 

}avq 

SR 

t{s} 

/best 

favg 

SR 

t{s) 

MDPIl 

3000 

189.048965 

189.048965 

20/20 

88.36 

189.048965 

189.048965 

20/20 

120.05 

MDPI2 

3000 

187.387292 

187.387292 

20/20 

60.71 

187.387292 

187.387292 

20/20 

101.07 

MDPI3 

3000 

185.666806 

185.655084 

13/20 

352.85 

185.666806 

185.651588 

10/20 

526.05 

MDPI4 

3000 

186.163727 

186.153631 

16/20 

300.37 

186.163727 

186.163727 

20/20 

136.64 

MDPI5 

3000 

187.545515 

187.545515 

20/20 

61.29 

187.545515 

187.545515 

20/20 

133.70 

MDPI6 

3000 

189.431257 

189.431257 

20/20 

51.99 

189.431257 

189.431257 

20/20 

35.59 

MDPI7 

3000 

188.242583 

188.242583 

20/20 

86.57 

188.242583 

188.242583 

20/20 

137.56 

MDPI8 

3000 

186.796814 

186.796814 

20/20 

48.04 

186.796814 

186.796814 

20/20 

66.76 

MDPI9 

3000 

188.231264 

188.231264 

20/20 

151.78 

188.231264 

188.231264 

20/20 

101.11 

MDPIIO 

3000 

185.682511 

185.623778 

10/20 

228.72 

185.682511 

185.672371 

18/20 

352.70 

MDPIIl 

3000 

252.320433 

252.320433 

20/20 

59.70 

252.320433 

252.320433 

20/20 

72.49 

MDPII2 

3000 

250.062137 

250.062137 

20/20 

220.10 

250.062137 

250.054744 

7/20 

513.80 

MDPII3 

3000 

251.906270 

251.906270 

20/20 

146.32 

251.906270 

251.906270 

20/20 

127.32 

MDPII4 

3000 

253.941007 

253.940596 

19/20 

370.76 

253.941007 

253.939680 

18/20 

352.64 

MDPII5 

3000 

253.260423 

253.260350 

17/20 

374.00 

253.260423 

253.260164 

14/20 

349.19 

MDPII6 

3000 

250.677750 

250.677750 

20/20 

55.35 

250.677750 

250.677750 

20/20 

69.78 

MDPII7 

3000 

251.134413 

251.134413 

20/20 

74.72 

251.134413 

251.134413 

20/20 

97.74 

MDPII8 

3000 

252.999648 

252.999648 

20/20 

79.82 

252.999648 

252.999648 

20/20 

115.84 

MDPII9 

3000 

252.425770 

252.425770 

20/20 

90.27 

252.425770 

252.425770 

20/20 

106.79 

MDPIIIO 3000 

252.396590 

252.396590 

20/20 

13.18 

252.396590 

252.396590 

20/20 

16.06 

T^Better 


0 

4 

4 

14 

0 

2 

2 

6 

T^Equal 


20 

14 

14 

0 

20 

14 

14 

0 

9 ^ Worse 


0 

2 

2 

6 

0 

4 

4 

14 

p-value 


1.0 

4.142e-l 







MDPIl 

5000 

240.162535 

240.102875 

7/20 

312.13 

240.141212 

240.021201 

2/20 

163.67 

MDPI2 

5000 

241.827401 

241.792978 

6/20 

1244.36 

241.817543 

241.753546 

2/20 

734.33 

MDPI3 

5000 

240.890819 

240.888162 

19/20 

810.48 

240.890819 

240.825167 

4/20 

531.94 

MDPI4 

5000 

240.997186 

240.976789 

6/20 

653.64 

240.973489 

240.915459 

3/20 

560.43 

MDPI5 

5000 

242.480129 

242.475885 

19/20 

735.16 

242.480129 

242.430474 

5/20 

515.62 

MDPI6 

5000 

240.322850 

240.306326 

8/20 

976.02 

240.328684 

240.266264 

5/20 

290.72 

MDPI7 

5000 

242.814943 

242.774982 

5/20 

259.50 

242.820139 

242.759895 

3/20 

819.01 

MDPI8 

5000 

241.194990 

241.161763 

8/20 

1148.60 

241.144781 

241.113453 

3/20 

721.59 

MDPI9 

5000 

239.760560 

239.667613 

4/20 

1219.71 

239.760560 

239.514958 

4/20 

360.05 

MDPIIO 

5000 

243.473734 

243.373015 

4/20 

457.28 

243.385487 

243.348149 

9/20 

939.59 

MDPIIl 

5000 

322.235897 

322.181291 

5/20 

1519.05 

322.223220 

322.131204 

3/20 

197.42 

MDPII2 

5000 

327.301910 

327.006342 

5/20 

1103.13 

327.301910 

327.075247 

5/20 

18.52 

MDPII3 

5000 

324.813456 

324.801590 

10/20 

955.81 

324.810826 

324.790223 

4/20 

27.82 

MDPII4 

5000 

322.237586 

322.197276 

5/20 

664.10 

322.212289 

322.126605 

4/20 

338.02 

MDPII5 

5000 

322.491211 

322.380726 

7/20 

1014.90 

322.420806 

322.301249 

5/20 

105.12 

MDPII6 

5000 

322.950488 

322.703887 

4/20 

352.88 

322.950488 

322.615227 

5/20 

212.55 

MDPII7 

5000 

322.850438 

322.793125 

10/20 

714.31 

322.850438 

322.778396 

8/20 

756.14 

MDPII8 

5000 

323.112120 

323.053268 

11/20 

879.48 

323.033840 

322.873156 

5/20 

161.94 

MDPIIO 

5000 

323.543775 

323.339842 

7/20 

569.73 

323.522709 

323.278556 

3/20 

148.94 

MDPIIIO 5000 

324.519908 

324.414458 

15/20 

752.95 

324.519908 

324.294790 

10/20 

753.14 

T^tBetter 


11 

19 



2 

1 



:;)tEqual 


7 

0 



7 

0 



# Worse 


2 

1 



11 

19 



p-value 


1.26e - 2 

5.699e - 5 








algorithm performs slightly better than the MTS algorithm, but the differences 
is small. However, for the 20 larger instances with n = 5000, the MAMMDP 
algorithm significantly outperforms the the MTS algorithm. First, compared 
with the MTS algorithm, the MAMMDP algorithm obtains better and worse 
results in terms of the best objective value on 11 and 2 instances respectively. 
Second, in terms of average objective value, the MAMMDP algorithm yields 
better results on 19 out of 20 instances. In addition, from the Friedman test, 
one observes that the obtained p-values are 1.26e — 2 (<0.05) and 5.699e — 5 
(<0.05) respectively for the best and average objective values, implying there 
exists a significant difference between these two methods. These outcomes 
demonstrate that the memetic framework is particularly useful to solve large 
and difficult instances. 
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5 Conclusions 


In this paper, we propose the first population-based memetic algorithm (MAM- 
MDP) for solving the NP-hard max-mean dispersion problem (MaxMeanDP). 
MAMMDP integrates an effective tabu search procedure and a random crossover 
operator while adopting an original scheme for parent selection. The compu¬ 
tational results on a large number of 200 benchmark instances show that the 
proposed algorithm is very competitive compared with the state-of-the-art al¬ 
gorithms in the literature. Specihcally, it improves or matches the previous 
best known results for all tested instances with n < 1000 with an average 
computing time of less than 14 seconds and a success rate of 100%, with only 
one exception. In particular, we found new and improved best results for 59 
out of the 60 most challenging instances. We also show computational results 
on 40 large instances with 3000 or 5000 elements which can serve as reference 
lower bounds for evaluating new MaxMeanDP algorithms. 

The investigations of several important ingredients confirm that both the un¬ 
derlying tabu search procedure and the crossover operator of the proposed 
algorithm contribute to the high performance of the proposed algorithm. It is 
shown that the population-based memetic framework is particularly suitable 
to solve large and difficult problem instances. 

The proposed algorithm could be adapted to the weighted version of the max- 
mean dispersion problem with several small modihcations. Some ideas of the 
proposed algorithm could be applied to other other binary optimization prob¬ 
lems (including some dispersion problems) where no constraint is imposed on 
the number of variables taking the value of one. 
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Table 6 

Appendix: Computational results of the proposed MAMMDP algorithm on the 
set of 100 small instances with n < 150. MAMMDP is run 20 times to solver each 
instance, each run being limited to 10 seconds. The previous best known results 
ifpre) are from: http://www.optsicom.es/edp/ 

Memetic Algorithm 


Instance 

n 

fvre 

/best 



t{s) 

MDPIl 20 

20 

13.88 

13.880000 

13.880000 

20/20 

0.02 

MDPI2 20 

20 

13.608 

13.608000 

13.608000 

20/20 

0.02 

MDPI3 20 

20 

11.7957 

11.795714 

11.795714 

20/20 

0.02 

MDPI4 20 

20 

17.54 

17.540000 

17.540000 

20/20 

0.02 

MDPI5 20 

20 

16.0063 

16.006250 

16.006250 

20/20 

0.02 

MDPI6 20 

20 

14.6064 

14.606364 

14.606364 

20/20 

0.02 

MDPI7 20 

20 

14.8822 

14.882222 

14.882222 

20/20 

0.02 

MDPI8 20 

20 

14.4614 

14.461429 

14.461429 

20/20 

0.02 

MDPI9 20 

20 

14.035 

14.035000 

14.035000 

20/20 

0.02 

MDPIIO 20 

20 

13.4433 

13.443333 

13.443333 

20/20 

0.02 

MDPIIl 20 

20 

18.855 

18.855000 

18.855000 

20/20 

0.02 

MDPII2 20 

20 

17.83 

17.830000 

17.830000 

20/20 

0.02 

MDPII3 20 

20 

18.11 

18.110000 

18.110000 

20/20 

0.02 

MDPII4 20 

20 

17.842 

17.842000 

17.842000 

20/20 

0.02 

MDPII5 20 

20 

16.344 

16.344000 

16.344000 

20/20 

0.02 

MDPII6 20 

20 

17.61 

17.610000 

17.610000 

20/20 

0.02 

MDPII7 20 

20 

18.9383 

18.938333 

18.938333 

20/20 

0.02 

MDPII8 20 

20 

21.88 

21.880000 

21.880000 

20/20 

0.02 

MDPII9 20 

20 

19.785 

19.785000 

19.785000 

20/20 

0.02 

MDPIIIO 20 

20 

22.599 

22.599000 

22.599000 

20/20 

0.02 

MDPIl 25 

25 

17.2708 

17.270833 

17.270833 

20/20 

0.03 

MDPI2 25 

25 

15.1214 

15.121429 

15.121429 

20/20 

0.02 

MDPI3 25 

25 

14.1817 

14.181667 

14.181667 

20/20 

0.03 

MDPI4 25 

25 

19.8567 

19.856667 

19.856667 

20/20 

0.03 

MDPI5 25 

25 

17.5371 

17.537143 

17.537143 

20/20 

0.02 

MDPI6 25 

25 

17.9667 

17.966667 

17.966667 

20/20 

0.03 

MDPI7 25 

25 

16.207 

16.207000 

16.207000 

20/20 

0.03 

MDPI8 25 

25 

18.1367 

18.136667 

18.136667 

20/20 

0.03 

MDPI9 25 

25 

17.4778 

17.477778 

17.477778 

20/20 

0.03 

MDPIIO 25 

25 

19.4592 

19.459167 

19.459167 

20/20 

0.03 

MDPIIl 25 

25 

21.81 

21.810000 

21.810000 

20/20 

0.03 

MDPII2 25 

25 

22.185 

22.185000 

22.185000 

20/20 

0.03 

MDPII3 25 

25 

23.5644 

23.564444 

23.564444 

20/20 

0.03 

MDPII4 25 

25 

19.74 

19.740000 

19.740000 

20/20 

0.03 

MDPII6 25 

25 

20.1744 

20.174444 

20.174444 

20/20 

0.03 

MDPII7 25 

25 

19.947 

19.947000 

19.947000 

20/20 

0.03 

MDPII8 25 

25 

23.921 

23.921000 

23.921000 

20/20 

0.03 

MDPII9 25 

25 

25.016 

25.016000 

25.016000 

20/20 

0.03 

MDPIIIO 25 

25 

23.575 

23.575000 

23.575000 

20/20 

0.03 

MDPIl 30 

30 

19.861 

19.861250 

19.861250 

20/20 

0.03 

MDPI2 30 

30 

18.813 

18.813333 

18.813333 

20/20 

0.03 

MDPI3 30 

30 

15.249 

15.248889 

15.248889 

20/20 

0.03 

MDPI4 30 

30 

22.717 

22.717333 

22.717333 

20/20 

0.03 

MDPI5 30 

30 

17.237 

17.236667 

17.236667 

20/20 

0.03 

MDPI6 30 

30 

18.375 

18.375455 

18.375455 

20/20 

0.03 

MDPI7 30 

30 

15.293 

15.292500 

15.292500 

20/20 

0.03 

MDPI8 30 

30 

19.247 

19.247273 

19.247273 

20/20 

0.04 

MDPI9 30 

30 

22.004 

22.004286 

22.004286 

20/20 

0.03 

MDPIIO 30 

30 

18.698 

18.698462 

18.698462 

20/20 

0.03 

MDPIIl 30 

30 

22.2721 

22.272143 

22.272143 

20/20 

0.03 

MDPII2 30 

30 

26.9138 

26.913846 

26.913846 

20/20 

0.03 

MDPII3 30 

30 

21.8973 

21.897273 

21.897273 

20/20 

0.04 

MDPII4 30 

30 

20.5375 

20.537500 

20.537500 

20/20 

0.03 

MDPII5 30 

30 

22.79 

22.790000 

22.790000 

20/20 

0.03 

MDPII6 30 

30 

20.351 

20.351000 

20.351000 

20/20 

0.03 

MDPII7 30 

30 

27.655 

27.655000 

27.655000 

20/20 

0.03 

MDPII8 30 

30 

26.8842 

26.884167 

26.884167 

20/20 

0.03 

MDPII9 30 

30 

24.1767 

24.176667 

24.176667 

20/20 

0.03 

MDPIIIO 30 

30 

24.8 

24.800000 

24.800000 

20/20 

0.03 

MDPIl 35 

35 

19.1833 

19.183333 

19.183333 

20/20 

0.04 

MDPI2 35 

35 

17.168 

17.168000 

17.168000 

20/20 

0.04 

MDPI3 35 

35 

17.0746 

17.074615 

17.074615 

20/20 

0.04 

MDPI4 35 

35 

23.35 

23.350000 

23.350000 

20/20 

0.04 

MDPI5 35 

35 

19.0177 

19.017692 

19.017692 

20/20 

0.04 

MDPI6 35 

35 

19.445 

19.445000 

19.445000 

20/20 

0.04 

MDPI7 35 

35 

19.4971 

19.497143 

19.497143 

20/20 

0.04 

MDPI8 35 

35 

21.2307 

21.230667 

21.230667 

20/20 

0.04 

MDPI9 35 

35 

20.98 

20.980000 

20.980000 

20/20 

0.04 

MDPIIO 35 

35 

16.9378 

16.937778 

16.937778 

20/20 

0.04 

MDPIIl 35 

35 

25.968 

25.967500 

25.967500 

20/20 

0.04 

MDPII2 35 

35 

26.135 

26.135455 

26.135455 

20/20 

0.04 

MDPII3 35 

35 

24.159 

24.159231 

24.159231 

20/20 

0.04 

MDPII4 35 

35 

24.415 

24.415000 

24.415000 

20/20 

0.04 

MDPII5 35 

35 

23.857 

23.857500 

23.857500 

20/20 

0.04 

MDPII6 35 

35 

24.673 

24.673077 

24.673077 

20/20 

0.04 

MDPII7 35 

35 

29.394 

29.393846 

29.393846 

20/20 

0.04 

MDPII8 35 

35 

25.297 

25.217273 

25.217273 

20/20 

0.04 

MDPII9 35 

35 

27.435 

27.435000 

27.435000 

20/20 

0.04 

MDPIIIO 35 

35 

25.712 

25.712500 

25.712500 

20/20 

0.04 

MDPIl 150 

150 

45.92 

45.920192 

45.920192 

20/20 

0.17 

MDPI2 150 

150 

43.39 

43.392381 

43.392381 

20/20 

0.19 

MDPI3 150 

150 

40.05 

40.046304 

40.046304 

20/20 

0.17 

MDPI4 150 

150 

44.04 

44.044138 

44.044138 

20/20 

0.17 

MDPI5 150 

150 

42.48 

42.479388 

42.479388 

20/20 

0.17 

MDPI6 150 

150 

43.72 

43.722955 

43.722955 

20/20 

0.17 

MDPI7 150 

150 

46.08 

46.077308 

46.077308 

20/20 

0.16 

MDPI8 150 

150 

42.45 

42.451346 

42.451346 

20/20 

0.18 

MDPI9 150 

150 

42.48 

42.479767 

42.479767 

20/20 

0.18 

MDPIIO 150 

150 

41.8 

41.797805 

41.797805 

20/20 

0.16 

MDPIIl 150 

150 

57.48 

57.484000 

57.484000 

20/20 

0.16 

MDPII2 150 

150 

57.82 

57 . 4^^52 

57.820652 

20/20 

0.17 

MDPII3 150 

150 

58.42 

58.421818 

58.421818 

20/20 

0.17 

MDPII4 150 

150 

57.38 

57.381064 

57.381064 

20/20 

0.17 

MDPII5 150 

150 

54.23 

54.228571 

54.228571 

20/20 

0.17 

MDPII6 150 

150 

56.44 

56.442653 

56.442653 

20/20 

0.16 

MDPII7 150 

150 

58.89 

58.889167 

58.889167 

20/20 

0.17 






